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Abstract 

Recently, a method was proposed and tested to find saddle points of the action in simula- 
tions of non-abelian lattice gauge theory. The idea, called 'extremization', is to minimize 
J{SS/6A^y. The method was implemented in an explicitly gauge variant way, however, 
and gauge dependence showed up in the results. 

Here we show how extremization can be formulated in a way that preserves gauge 
invariance on the lattice. The method applies to any gauge group and any lattice action. 
The procedure is worked out in detail for the standard plaquette action with gauge groups 
U(l) and SU(N). 
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1 Introduction 



One of the approaches to understand the non-perturbative nature of non-abehan gauge 
theories has been to look for topological objects like monopoles and instantons in lattice 
Monte Carlo simulations. They are local minima of the action and are thought to give 
important contributions to the path integral. An important question is whether the 
presence of such objects is related to the string tension. The problem one faces here is 
that short distance fluctuations obscure the presence of such topological objects in lattice 
gauge field configurations. Monopoles and instantons are basically semiclassical objects, 
characterized by their long distance structure. 

Several techniques have been developed to extract the relevant long distance infor- 
mation from a lattice gauge field configuration. One of them is the cooling method 
which minimizes the action density of a particular configuration locally in an iterative way, 
directing it towards a local minimum in configuration space. Apart from local minima, 
however, saddle points are also solutions to the field equations. The sphaleron is a popular 
example. It is therefore useful to consider saddle points in lattice field configurations as 
well. 

Recently, Duncan and Mawhinney have discussed a method to find saddle point so- 
lutions of Yang- Mills theory (or other (lattice) field theories). They applied their method, 
which they called 'extremization', to three-dimensional SU(2) gauge theory and found that 
configurations after extremization showed lumps in the action density that were highly 
correlated to lumps found after cooling the same configurations. It was found that the 
string tension survives moderate extremization but ultimately disappears. 

The basic idea of extremization is to minimize the squared gradient of the action, 
referred to as 'extremization action', 
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(1) 

instead of lowering the action itself as in the cooling method. In ref. |@] this was imple- 
mented for the plaquette action of lattice gauge theory. A conjugate gradient method 
in terms of the gauge fields ^^(x) was used to minimize S and Fourier acceleration was 
used to speed up the minimization process. The authors stated that S" is a gauge variant 
function, so that lattices differing by a gauge transformation might extremize differently, 
which is what they observed in their gauge variant minimization procedure. 
However, if we take for S the (continuum) Yang- Mills action 

= ' (2) 

we see that the functional derivative in eq. (P gives just the gauge covariant equation of 
motion for A^(x), so that the corresponding extremization action S is the gauge invariant 
quantity 

^cont ^ f p^^a^)2 ^ (3) 

disregarding the factor 1/ g'^ in front of 5™°*. This is not surprising, for the extrema of 
the action, determined by 5, should be gauge invariant. Since the continuum quantity 
eq. (^ is gauge invariant, it must be possible to define the lattice analogue in a gauge 
invariant way as well. It is the purpose of this paper to show how this can be achieved. 
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In section 2 the gauge invariant extremization action is defined. It is shown that it has 
the expected behavior in the limit of lattice spacing a — >■ and its structure is visualized 
in terms of closed Wilson loops. Section 3 contains a discussion of the generality of the 
method and a few comments about its implementation. Finally, we briefly discuss physical 
applications. 



2 Gauge invariant extremization action 

Consider the standard plaquette action for SU(N) lattice gauge theory, 



S = /?E(l-^(Trf/n + Tr[/t) 

(/3 = 2N/g') 
= V (Tr f/n + Tr U^) + Const. 



(4) 



(5) 



We will forget about the irrelevant constant. The degrees of freedom in this theory are the 
link variables Un{x) which are matrices in a unitary representation of the gauge group. 
Therefore it is appropriate to consider the functional derivatives of the action with 
respect to U^{x). Since the link variable f/^(x) is related to a corresponding continuum 
gauge field A^{x) by the path ordered integral 



U = U^ix) 



—ia / An(x + tajl)dt 
Jo 



X 



X + a(x 



which can be expanded in a to give 

Uf,{x) = 1 - iaA^{x) + 0{a?) 



(6) 



(7) 



one expects this derivative to become equivalent to the derivative with respect to A^{x) 
in the limit a — > 0. In this way, we shall define a manifestly gauge invariant extremization 
action on the lattice and show that it reduces to eq. @ to leading order in the lattice 
spacing. 

We start by rewriting the plaquette action in a form that exhibits the dependence 
on a particular link U more clearly. 



S = -^ETr iUFu + U^Fl 



^ links 

u 



(8) 



(from now on we will omit factors of l/g^ 
'staples' for the link U, 



+ terms independent of U 

Here Fu is the sum of the 7 = 2{d — 1) 



/ 



u ■ 



E 



+ 



X 



(9) 



X 



3 



which upon contraction with U gives the part of the plaquette action depending on U. 7 
can be considered as the coordination number for the interactions. A staple is defined as 



f/^(a; + afi) UUx + az>) Ul{x) 



(10) 



X 



Fjj is sometimes called the 'force' (hence the notation F) for the link U . Note that it is 
not an element of the SU(N) in general. In the special case of gauge group SU(2) it can 
always be written as a constant times an SU(2) matrix. 

A variation 5U in the link U = Ufj,{x) gives rise to a change in action 6S equal to 



'6S 
6U 



5S = -Tr 

= -Tt [{UF - {UFy)6UU^} , 

where we have used that 

= 6{U^U) = {6U^)U + U^U , 
and it is understood that F = Fjj. We can visualize {UF — {UF)'^) as 



(11) 
(12) 

(13) 



E 



X 



X 



\ 



(14) 



where the sum runs over the 7 = 2((i — 1) staple directions. 

We now define the gauge invariant lattice analogue of the extremization action S by 



S = E^Tr {(f/F-(t/F)t)t(f/F-(f/F)t)} 



links 

u 



-Tr {2F^F - FUFU - {FUFU)^} . 

links 

u 



(15) 
(16) 



This is a sum over closed Wilson lines of lengths 6 and 8, arising from contractions of the 
diagrams in eq. ([l4|) . Obviously, S is gauge invariant. 

Before giving a more detailed survey of the diagrams we shall show that S reduces 
to its continuum counterpart (^ in the limit a — 0. For simplicity, we shall do so for 
the abelian case. The non-abelian generalization will follow easily. For completeness we 
mention that for the gauge group U(l) the plaquette action is 



U{1) 



and the appropriate definition of S is 



2r 



E(° 



(17) 



SuH) = \ EiUF - (UFmUF - (UFY) . 

links 

u 
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Consider the expression for UF—{UF)* in eq. (0) with U = U^{x). For each v ^ [i the 
sum contains a contribution in the positive and one in the negative //-direction. The term 
from the positive i/-direction is calculated as follows. First one computes the expansion 
of the plaquette around its centre, 



X 



■ exp [—ia? F^y{x + \ajl + |az>) + 0{a^)\ 



(19) 



so that 



X 



X 



V 



-2i sin [a?F^y{x 



Oia')] . (20) 



If the expansion in eq. (O) were not done around the centre of the plaquette, 0{a^) terms 
would be present as well. The contribution from the negative zz-direction is 



+ 2i sin [a'^F^^{x 



(21) 



and the sum of eqs. (^) and (^T]) can be expanded around x + |a/i and subsequently 
summed over v to give for UF — (UF)* in eq. (|I4D : 



/ 



E 



X 



X 



\ 



-2m=^ ^ a,F^,(x + ia/i) + C»(c 



(22) 



Terms of 0{a^) cancel out. Finally, this result is inserted in eq. ([I8|), giving 



5, 



C/(l) 



a'(9.F^.)' + 0(a' 



(23) 



which is the abelian version of eq. (^ to leading order in a, as promised. The extra factor 
of here is explained by eq. (|^) and serves to make S dimensionless. This result can be 
generalized to the non-abelian case immediately, by replacing the derivatives by covariant 
ones and adjusting the coefficients appropriately. 

Now we turn to an overview of the diagrams occurring in S (p!^-[T6|). Each link U in 
the summation eq. (p!6D yields (27)^ closed Wilson loops, with the appropriate signs (recall 
that 7 = 2{d — 1)). In 2, 3 and 4 dimensions, 47^ = 16, 64, 144 respectively. These huge 
numbers include a number of trivial and double-counted diagrams as well as hermitian 
conjugates, however. 

The first term in eq. ([T6| ) contributes 7^ Wilson loops of length 6 in lattice units, each 
counted twice. These are the diagrams 'hinging' around the link Un{x), but independent 
of it, together with their adjoints, see fig. la. Of these, 7 give the unit matrix, another 
7 are the planar 2x1 loops, and the remaining 7(7 — 2) do not lie in a plane. Note 
that the symmetric length 6 Wilson loops winding around a cube as depicted in fig. lb 
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Figure 1: 

a. Length 6 Wilson loop contributing to S. 

b. Non-contributing length 6 diagram. 

c. Contributing diagram of length 8. 

d. 'Double plaquette'. 

are excluded. The other terms in eq. (|16]) give the length 8 Wilson loops containing the 
link U twice (fig. Ic) and their conjugates. These 27^ diagrams include the 7 'double 
plaquettes' shown in fig. Id with their conjugates. In these double plaquette diagrams, 
the loop winds around a plaquette twice before closing on itself. 

Note how nicely this method provides one with a lattice action reducing to the contin- 
uum derivative action eq. @ in lowest order in the lattice spacing a. It is also interesting 
to observe that, apart from being gauge invariant, eq. (|^) preserves Lorentz-invariance. 
This implies that S cannot be used as an improvement term for the plaquette action. 
The 0{a^) term in the expansion of the plaquette action which is canceled in tree-level 
improved actions ^ by adding 2x1 rectangular Wilson loops, are of the form 

J2D,F,^D,F,^ (24) 

fJ.,U 

which does not equal Y^^i^Jlv ^vFu^)"^ ■ 



3 Further remarks 

It should be clear that the applicability of the procedure described in this paper is not 
restricted to the standard plaquette action. For any lattice action S*, written as a sum 
over closed Wilson lines, the quantity 



S 



-VTr 
2 ^ 



5S 



6U 



(25) 



provides a gauge invariant extremization action. The calculational recipe of 5* can be 
visualized as follows. Cut the link U out of the fZ-dependent closed loops occurring in S, 
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and paste the resulting open-ended lines together (with the correct signs) to form closed 
loops of S. 

Furthermore, the extremization action eq. (^) can be defined for any gauge group, 
including discrete groups such as Z^, the link variables being group elements. Of course, 
it is not possible to define a continuum limit of the extremization action (nor the action 
itself) if the gauge group is discrete. This does not limit the use of the extremization 
action for lattice gauge theories with discrete groups, however. 

We should also make some remarks about implementation of the suggested extremiza- 
tion procedure. First of all, in four dimensions the number of diagrams contributing to S 
for each link is large so that calculation and minimization of S can be computationally 
quite expensive. In two and three dimensions this problem is less severe and the method 
may find its applications. Furthermore, minimization of the extremization action S is not 
straightforward. The dependence of on a particular link U is of the non-linear form 

Tr [MUMU + NU + hermitian conjugate] , (26) 

where both M and are sums over SU{N) matrices (or U(l) numbers). An iterative 
procedure is probably required to carry out the local minimization step. Finally, one may 
need special algorithms to speed up the long distance modes, as in ref. f^. 

We end with a few remarks concerning the application of extremization to physical 
problems, in the light of the work of Duncan and Mawhinney 0. It would be interesting 
to see if their results are confirmed when the gauge invariant procedure proposed here is 
used to extremize configurations. In any case, gauge dependence in the results should be 
absent when the gauge invariant procedure is used. 

In the simulation of three-dimensional SU(2) gauge theory in ref. [0 it was found that 
during extremization the action itself decreased as well. At first sight, this is somewhat 
surprising, as it is hard to see why a quantity should decrease when its derivatives are 
minimized. One wonders whether the fact that the functional derivatives in the gauge 
variant extremization procedure are not exact lattice derivatives has any influence here. 
Similarly, it is surprising that a strong correlation was noticed between lumps in the 
action density after cooling a configuration and after extremizing the same configuration. 
Finally, it was observed that the string tension disappeared after continued extremization. 
These are all interesting problems that merit further investigation. 
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